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Abstract 

We compute the condensate in QCD with a single quark flavor using numerical simulations with 
the overlap formulation of lattice fermions. The condensate is extracted by fitting the distribution of 
low lying eigenvalues of the Dirac operator in sectors of fixed topological charge to the predictions of 
Random Matrix Theory. Our results are in excellent agreement with estimates from the orientifold 
large- N c expansion. 
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I. INTRODUCTION 



Very few analytic techniques are available to study nonperturbative properties of QCD. 
Of these, the most prominent are large-iV c expansions. Recently, Armoni, Shifman, and 
Veneziano suggested a new large- N c expansion with some remarkable features. In 

contrast to the 't Hooft large- iV c limit [5( (N c — » oo, g 2 N c and Nf fixed, with quarks in the 
fundamental representation of SU(N C )), quarks are placed in the two- index antisymmetric 
representation of SU(N C ). Now in the N c — ► oo, g 2 N c and Nf fixed limit of QCD, quark 
effects are not decoupled, because there are as many quark degrees of freedom as gluonic 
ones, 0{N 2 ) in either case. In Ref. j3| the authors have argued that a bosonic sector of 
M = 1 super- Yang-Mills theory is equivalent to this theory in the large- N c limit. (The proof 
of this connection has recently been extended to lattice regularized theories by Patella 
The large-iV c QCD-like theory is called "orientifold QCD." 

For N c = 3, orientifold QCD is equivalent to QCD with a single quark flavor in the 
fundamental representation of SU (3). This equivalence even extends to the first and second 
terms in the (3 function and in the lowest order anomalous dimension for the running quark 
mass (or quark condensate). This means that nonperturbative quantities computed in super- 
Yang-Mills theory can be related to corresponding ones in one-flavor QCD, up to 1/N C effects. 

In a recent paper [j], Armoni, Shifman and Veneziano estimate the quark condensate in 
one-flavor QCD from the value of the gluino condensate in SUSY Yang-Mills. They find 
(with our sign conventions) 



in the MS scheme at \x = 2 GeV. The spread of values gives their estimate of 1/N C cor- 
rections. This corresponds to S 1 ' 3 of 240 to 300 MeV. They estimate values of S from the 
Gell-Mann, Oakes, Renner relation and from an interpolation of lattice data, which are in 
good agreement with Eq. ([I]). However, a direct lattice calculation of the condensate in 
Nf = 1 QCD would give a more reliable comparison. Such a calculation would be impor- 
tant to researchers studying the orientifold theory, for it would indicate the size of l/N c 
corrections to calculated quantities. This lattice calculation we now provide. 

In the literature, there are different but related quantities called the quark condensate. 
Often, the expression refers to (qq), which in a lattice simulation, as well as in the continuum, 
is a function of the quark mass and the volume in which the quark fields are defined. 

In this paper, however, we attempt to extract the low-energy constant S, i.e. a parameter 
of the low-energy effective theory. The case of one flavor is a bit special: the chiral symmetry 
is anomalous, there is a massive pseudoscalar, the r/, and no Goldstone bosons. The E which 
we are about to extract is therefore not an order parameter of spontaneous chiral symmetry 
breaking. However, there still exists a well defined low-energy description of Nf = 1 QCD. 
It has been worked out by Leutwyler and Smilga 0] to which we refer the reader for details. 
They show that up to terms of order m 2 V the partition function is 



£ = {0.014, 0.021, 0.028} GeV 3 



(1) 



Z = exp {Y*V~Re(me 



)} 



(2) 
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with 9 the vacuum angle. E is the infinite volume zero quark mass limit of — (qq) at 9 = 0. 
The lattice calculation has four parts: 

First, we must describe how to perform simulations with one dynamical quark flavor. 
These require the use of overlap fermions, for which algorithms have only recently been 
developed |. Regularities in the spectrum of overlap fermions allow us to use a version of 
the Hybrid Monte Carlo algorithm with chiral pseudofermions to give us an exact algorithm 
for any number of flavors [9, The use of a lattice action with exact chiral symmetry 
means that observables (like the anomaly) are not contaminated by explicit chiral symmetry 
breaking effects from the discretization. 

Since direct measurements of (qq) are influenced by finite volume and non-zero quark 
mass, we use another quantity whose dependence on the condensate is known. For this we 
choose the low eigenvalues of the Dirac operator, measured in sectors of fixed topology v in 
a simulation volume V. The probabilit y d istribution of individual eigenvalues A n is given 



by Random Matrix Theory (RMT) ll|, ll2|, Il3( as a function of the dimensionless quantity 
\ n T,V, which depends parametrically on the combination m q IlV and, of course, Jdie number 



of flavors Nf. (Formulas for the special case of Nj = 1 can be found in Ref. [7|.) We use 
the specific method and predictions from Refs. [3, [3] , where all the relevant formulas are 
displayed. 

Third, we need a lattice spacing to convert the dimensionless lattice-regulated condensate 
to a dimensionful number. The calculation of Ref. JsJ] uses a A parameter to set the scale. 
Lattice calculations have not done that for many years, since typically it is hard to extract a 
A parameter from simulation data and because most observed quantities in fact are not very 
sensitive to it. Instead, it is customary to set the scale with some spectral quantity, which 
can be computed directly in a simulation. Possibilities include the masses of various mesons 
or the string tension. We follow common lattice practice and obtain the lattice spacing 
through the Sommer parameter tq fig , which is defined through the force, 



r2 dV(r) 



Or 



r=r 



1.65 . (3) 



We could alternatively use the string tension, but r is less noisy. 

The prediction of Eq. flTJ is a number in GeV. Since Nj = 1 QCD has no physical 
realization, its overall scale is unknown. We will convert our dimensionless number into a 
dimensionful one using the real world value of r$ = 0.5 fm, even though its ratio to other 
observables is almost certainly Nf— dependent. 

Finally, we need a matching factor, to convert the lattice-regulated condensate to its MS 
value. We do this using; the Regulation Independent seherue 0. 

In Sec. Hi] we describe all the ingredients of the lattice calculation. Our results are 
summarized in Sec. inland we conclude in Sec. HVl 

The only previous calculation of the condensate for Nf = 1 we are aware of is that of 
Ref. fl-jj: the authors only quote Sa 3 . The calculation was done deep in the strong coupling 
limit with a single flavor of staggered fermions. Flavor symmetry was so badly broken that 
the four tastes act as a single physical flavor. 
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II. LATTICE METHODOLOGY 



A. Lattice action and simulation parameters 



Our simulations are performed with overlap fermions lg, 119J . This discretization of the 
Dirac operator preserves exact chiral symmetry at nonzero lattice spacing via the Ginsparg- 
Wilson relation [20|]. The massless operator is 



D = D r 



m 



0)=R [l+ 7B£ (/ l (-i2 ))] 



(4) 



with e(h) = h/yh? the sign function of the Hermitian kernel operator h = 75c/ which is taken 
at negative mass Rq. It has a spectrum consisting of chiral modes with real eigenvalues 
(at A = and 2Rq) and nonchiral modes which are paired complex conjugates, A and 
A*. The squared Hermitian overlap operator H 2 = (75.D) 2 = D^D commutes with 75 and 
therefore its eigenvectors can be chosen with definite chirality. Because of the Ginsparg- 
Wilson relation, to each eigenvalue |A| 2 correspond two eigenmodes of opposite chirality, 
and the eigenvectors of D ov {m = 0) with complex conjugate eigenvalues are superpositions 
of these two eigenvectors. It is convenient to define the chiral projections (P± = |(1 ± 75)) 
so that the massive squared Hermitian overlap operator, with the usual convention for the 
mass terms, is 



Hl(m) = P±H 2 {m)P± 



2{R 2 



m 



^)P ± (l±e(h))P ± +m 2 P ± . 



(5) 



Since the spectrum is doubled, the spectrum of one chiral sector of H 2 is equal to the 
spectrum of a single quark flavor, apart from the real modes. Their contribution can be 
included exactly [9]. Choosing to work in the chiral sector a which has no zero modes, the 
partition function for a single flavor of quark is then 



[d(j>l][d4a\ exp(-4 H a (m) 2 <p a - \Q\ log(m/(2i? ))). 



(6) 



where Q is the topological charge as defined by the number of zero modes of negative and 
positive chirality, Q = n_ — n + . This systemjs then simulated by the Hybrid Monte Carlo 
algorithm 



21 



with the modification of Ref. [8j for topological boundaries 



As we have described in Ref. lOj, the pseudofermion fields are initialized with a set of 
chiral Gaussian random vectors £ CT by CT = a/ H a (m) 2 £ a . In our simulations we have used 
two variations of this algorithm. For computing E, we have restricted the simulation to 
fixed sectors of Q. To find the string tension and matching factor, we have fixed the running 
chirality to be negative and restricted allowed topologies to be Q < 0. In the analysis of an 
ensemble generated with this algorithm, measurements on the Q = configurations need 
to be reweighted with a factor 1/2 compared to those from configurations with non-trivial 
topology. 

Our particular implementation of the Hybrid Monte Carlo algorithm has been previously 
discussed in Refs. H Q E 



It uses multiple pseudofermions and stepsizes [25 



26 
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We are using a planar kernel Dirac operator d with nearest and next-to-nearest ( "y^" ) 
interactions. To be precise, we parameterize the associated massless free action by 

S = J^O) [v{r) +i'j l iPn(r)}4>(x + r) (7) 

x,r 

with r connecting nearest neighbors (r = ±p; rj = r)i, p M = pi) and diagonal neighbors 
(r = ±/i ± i>, 1/ 7^ p; rj = r)2, Pn = p v = p%)- The constraint rj(r = 0) = t)q = —8771 — 24r? 2 
enforces masslessness on the spectrum, and —1 = 2pi + 12p 2 normalizes the action to 
—ipi'j fl d t j i ip in the naive continuum limit. To speed up the code we require that each of the 
couplings of a fermion to its neighbors is a projector, proportional to l±n-7- This is a familiar 
trick for Wilson action simulations. For nearest neighbors, a projector action corresponds to 
the constraint 771 = p\ (up to signs) and for the diagonal neighbors, r? 2 = v / 2p2- The action 
we use in the simulations presented in this paper uses p\ = —1/6 and p 2 = —1/18. We also 
add a clover term with the tree-level clover coefficient appropriate to this action of 1.278. 
We set the radius of the Ginsparg- Wilson circle Ro to 1.2. 

Our kernel operator d is constructed from gauge links to which three levels of isotropic 



stout blocking 27]] have been applied. The blocking parameter p is set to 0.15. The sign 



function is computed using the Zolotarev approximation with an exact treatment of the 



low-lying eigenmodes |A) of h(—R ). We use the Liischer-Weisz gauge action [28j with the 
tadpole improved coefficients of Ref. [29J]. Instead of determining the fourth root of the 
plaquette expectation value uq = {{U p i) /3) 1 / 4 self-consistently, we set it to 0.86 for all our 
runs as we did in our previous publications. We simulate on 10 4 lattices at one value of the 
gauge coupling ft — 7.7 which we chose to be roughly at a lattice spacing of 0.16 fm. Our 
bare sea quark mass is am q = 0.05. We collected approximately 600 trajectories of data in 
fixed Q=0 and 1 sectors, and analyzed lattices from every fifth trajectory. 

The use of three steps of stout blocking results in a considerable saving of computer time 
compared to the two steps we have used previously. This is seen by comparing the number 
of inner Conjugate Gradient steps needed to evaluate the Zolotarev approximation to the 
sign function. In these simulations it is about 20. With two steps, at essentially the same 
values of lattice spacing and quark mass, it is about 47. 

The range of smearing of the gauge fields for N steps stout blocking is [30j (x 2 ) ~ 2pa 2 N. 
This corresponds to ((x/a) 2 ) 1 / 2 ~ 0.77 for N = 2 and 0.95 for N = 3 at p = 0.15 such that 
also for the N = 3 case the fermion action remains reasonably local. 



B. The lattice-regulated condensate 



We computed the lowest four eigenvalues |A| 2 of the squared Dirac operator H 2 . They 
give us the eig envalues of the overlap operator, which lie on a circle. We apply the Mobius 



transform 31 



A 



A 



1 - A/(2i2o) 

to project the eigenvalues onto the imaginary axis. 



(8) 
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FIG. 1: Cumulative distributions of the eigenvalues where C(A) is the fraction of all gauge con- 
figurations where the n-th eigenvalue is smaller than |A|. We show the data for n = 1 and 2 for 
topological sector v = and 1. The dotted lines are the predictions from RMT with the value of 
Ea 3 obtained from fits to the data (from left to right: fits C, B, A from Table [J). 



The data set thus consists of a collection of eigenvalues, the distribution of which is 
predicted by RMT as presented in Refs. 0,0. It depends on one parameter E. 

To avoid binning the data, which can introduce a bias, we use the Kolmogorov-Smirnov 
test 3l|, |32| as a measure for the goodness of the fit. It compares the cumulative distribution 
function of the data C(x) to the theoretical prediction P(x) = f(x)dx. C(x) is the 
fraction of eigenvalues with a value smaller than x. 

The quantity of interest is the largest deviation of P and C: D = max x \P(x) — C(x)\. 
From this the confidence level is given by 

Qks((VN + 0.12 + 0.11/VN)D) (9) 

with 

oo 

Q KS (X) = 2^(-r 1 exp(-2j 2 A 2 ). (10) 

3=1 

In fits to a single eigenvalue distribution we maximize this quantity. When fitting to more 
than one eigenvalue, we maximize the product over the individual confidence levels. The 
errors on the fit parameter Sa 3 are determined by the bootstrap procedure. 

The results of our fits are presented in Table [land Figs. [TJ EJ From a fit to the lowest level 
in v = alone (fit A) we extract Sa 3 = 0.0087(4) with a confidence level(CL) of 0.87. This 
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0.0102(4) 
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0.59 




0.21 



TABLE I: The confidence level of the individual distributions from the various fits. The values for 
levels included in the fits are boxed. 




FIG. 2: Binned eigenvalue distributions and the RMT prediction with Sa 3 from the fit to the 
lowest level in each topological sector (fit B from the table). 

means that the shape of the extracted curve is in very good agreement with the theoretical 
prediction for the given statistics. The second level in v = has a CL of 0.02, which is 
low but still acceptable. A fit to the lowest level in \u\ = 1 (fit C) gives a 3 £a 3 = 0.0102(4) 
and a CL for the lowest two level of 0.59 and 0.21. From a combined fit to the lowest 
level in each sector (fit B) we get Sa 3 = 0.0096(3). Even with the extracted values of Sa 3 
seemingly apart, it is still possible that they are compatible with another, because there is 
no correlation between the v — and \u\ = 1 ensembles. 

Now that we have determined the optimal Sa 3 , the quality of the fits has to tell us 
whether the RMT description of the data is valid. A possible source for a deviation from 
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the RMT prediction is a too small volume. Since we simulate at finite volume, we only 
expect the lowest few eigenvalues to match the RMT curves. We observe that for the fits to 
one level alone, the fit to the second eigenvalue in the same topological sector still has an 
acceptable confidence level. The match between the £ in v = and \v\ — 1 does not seem 
convincing. However, one has to keep in mind that the lowest two eigenvalues are calculated 
on the same gauge configurations and therefore are correlated. On the other hand, there is 
no correlation between the extraction in the two topological sectors. This can give rise to an 
apparent better match between the distributions of the two levels in one topological sector 
as compared to the match between the lowest level in both sectors. Within our statistics, 
RMT seems to be applicable to the lowest eigenvalue in v = and \v\ — 1 and the second 
lowest eigenvalue in v = 0. Nevertheless, it is well possible that the volume is still too small. 

A main concern in any lattice calculation — and in particular ones using new actions and 
algorithms — is the auto-correlation between consecutive configurations. From our Markov 
chains, we saved every fifth configuration. To find out whether this separation is enough, we 
used the following technique: We split the data set into two sub-sets, one consisting of all the 
even numbered configurations and one of the odd numbered. A fit to the lowest eigenvalue 
in v = and \u\ = 1 on the even and the odd numbered sample gives a 3 X = 0.0099(5) and 
0.0094(4) respectively. The central values are almost one sigma apart suggesting a weak 
correlation. 

To make a more quantitative statement, we have computed a 3 £ on matched Bootstrap 
samples (where if configuration 2i is part of the "even" Bootstrap sample, so is 2i + 1 in the 
odd sample). Correlations between the two samples should show up as correlations between 
the £ on the corresponding bootstrap samples. We therefore measured the correlation matrix 
element (averaging over all bootstrap samples) 

C = ((a 3 E even — a 3 S e ven)(a 3 S dd — a 3 E dd)) (11) 

and found C/ (a 6 £ e ven£odd) = 0.0007(15) which is zero within errors. Thus, we do not 
find any effects of auto-correlation beyond the five trajectories by which we spaced our 
measurements. 

In a second test of autocorrelations, we computed the integrated autocorrelation time 
for individual eigenvalues from the data stream. We found integrated autocorrelation times 
averaging about 1.5 (in units of the collection time, 5 trajectories), with an uncertainty 
of about 0.7. This again suggests a weak correlation between successive measurements of 
eigenvalues. From this analysis we conclude that the value of the condensate in lattice units 
is a 3 S = 0.0096(3). 



C. The length scale 

As discussed above, an overall scale is obtained from the static quark potential. The latter 
is extracted from the effective masses of Wilson loops after one level of HYP smearing 33, 



34j, where the short-distance effects of the HYP smearing are corrected using a fit to the 
perturbative lattice artifacts. We measured the potential on 200 8 3 x 12 configurations at 
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FIG. 3: The static quark potential in lattice units. The filled symbols denote the potential after 
removing the artifacts introduced by the HYP smearing. 

am q = 0.05 with the result shown in Fig. [3l From the parameters of the fit we obtain 
ro/a = 3.37(10) for the Sommer parameter ([3]) and a\fa = 0.318(25) for the string tension. 
With r = 0.5 fm, this gives a lattice spacing of a = 0.15 fm. 



D. Matching lattice and continuum regularizations 

A matching factor is needed to convert the lattice calculation of the quark condensate to 
its corresponding MS value. To get the matching factor, we use the RT scheme introduced 



in Ref. [171], and we follow the procedure described in Ref. [35j, in which we calculated the 
matching factors for a quenched simulation. The RF scheme results in the chiral limit can 
be converted to MS values at fi = 2 GeV by using the ratios connecting the two schemes. 



The ratios were computed by continuum perturbation theory to three loops [36j, [37 . 

These simulations should not be restricted in topological sectors. To produce a matching 
factor, one needs simulations with a momentum scale short enough to be free from nonper- 
turbative effects and yet not so short as to be affected by discretization errors. Our data set 
consisted of 57 8 4 configurations from about 300 trajectories at am q = 0.05. 

The 8 4 lattice is periodic in space directions and antiperiodic in the time direction. There- 
fore the momentum values are 



[jk X} —k y} —k z ,-(2k t + l)y (12) 



We choose the values of k^ such that the momentum values lie as close as possible to 
the diagonal of the Brillouin zone. The maximum value of ap = 2.256 corresponds to 
kfj, = (2, 1, 1, 1). The quark propagators are cast from a point source and then projected to 
the desired momentum values. 
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FIG. 4: Zf v for the one flavor overlap simulation. 

TABLE II: Values of Z$ in the RF scheme. The inverse lattice spacing is 1.330 GeV from the 
Sommer parameter. Therefore /U = 2 GeV corresponds to a/i = 1.504. The value at this point is 
obtained from a linear interpolation from the two closest a/i points of the data. 

o/x 0.878 1.178 1.416 1.619 
Zf 0.76(6) 0.72(4) 0.73(3) 0.78(3) 

a/x 1.800 1.963 2.115 2.256 
Zf' 0.78(2) 0.772(15) 0.777(17) 0.780(11) 



Zg 1 ' for the scalar density is shown in Fig. HI The values of Zip' are listed in Table [TT1 
From our lattice spacing determined from the Sommer parameter, fx = 2 GeV corresponds 
to a fx = 1.504. The 2 GeV RF value is obtained from a linear interpolation from the two 
closest an points of the data. The result is Zf 1 ' (2 GeV) = 0.75(3). 

The conversion ratio from the RF to the MS scheme for the scalar and pseudoscalar 
densities, from [36|, |37|, needs an a s . In Landau gauge and to three loops, and for one flavor, 
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the ratio is given by 



Zf s Zf^_ 16 a, / 1375 _ 1526 \ / a*\ 2 

Z* 1 ' + 3 4vr + V6 3 J UJ 

/32149271 _ 2154896 _ 80^ 2960^ ,«.y 

V 2916 54 3 9 7 V4vr/ v s; ' v ; 

where Cn is the Riemann zeta function evaluated at n. 

To get numerical results of the above ratio, we use the coupling constant from the so- 



called "ay" scheme. As in the appendix of Ref. 23], from the one-loop expression relating 
the plaquette to the coupling 

\n^U p = -^a v {q*)W, (14) 

where W = 0.366 and q*a = 3.32 for the tree-level Liischer Weisz action, we obtain 
«y(3.32/a) = 0.173. Then aAj^ is calculated and a^ s (2 GeV) is determined by using 
& = 31/12tt and ft = 268/48vr 2 forone flavor QCD. We find af E (2_GeV) = 0.194. Substi- 
tuting af* into Eq. we get Z^/Zf = 1.147 and therefore Z^(2 GeV) = 0.86(3). 



III. RESULTS 



From Sec. [II] the continuum-regularized condensate is 



,3 w / r 0\3 



rZE(MS,fi= 2 GeV) = ZJfi,a) x Ea 3 x (— ) 

a 

= 0.86(3) x 0.0096(3) x (3.37(10)) 



3 



= 0.317(32). (15) 
Taking the real- world value for r = 0.5 fm, this is 

E(MS, fi = 2 GeV) = 0.0194(20)GeV 3 (16) 

or 

(E(MS, /x = 2 GeV)) 1/3 = 0.269(9)GeV. (17) 



IV. CONCLUSIONS 

Our result is in remarkable agreement with the central value of Eq. ([T]) . Thus at least for 
the condensate the 1 ± 1/N C estimates of the correction to the large N c result, the extreme 
values of Eq. (1T1), seem to be quite pessimistic. In view of this finding more predictions from 
large N c orientifold QCD might be interesting even without the knowledge of the subleading 
corrections. 

A summary of previous calculations of the condensate has recently been given by Mc- 



Neile 38j. Comparing the quenched determinations there, a three-flavor prediction by Mc- 



Neile, and our result, the condensate seems to be a quantity which is not very Nf dependent. 
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Other tests of orientifold equivalence will be difficult. To check the prediction of Ref. [l| 
that mir/rri^ = 1 + 0(1 /N c ) is nontrivial because it requires disconnected diagrams. A 
similar degeneracy of hybrids described in Ref. [isj will be hard because the sources for 
ordinary hybrids are noisy in QCD. 
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